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We study classical and quantum Heisenberg antiferromagnets with exchange anisotropy of XXZ-type and 
crystal field single-ion terms of quadratic and cubic form in a field. The magnets display a variety of phases, 
including the spin— flop (or, in the quantum case, spin-liquid) and biconical (corresponding, in the quantum 
lattice gas description, to supersolid) phases. Applying ground— state considerations, Monte Carlo and density 
matrix renormalization group methods, the impact of quantum effects and lattice dimension is analysed. 
Interesting critical and multicritical behaviour may occur at quantum and thermal phase transitions. 
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1. Introduction 

Uniaxially anisotropic Heisenberg antiferromagnets in a magnetic field along the easy axis have 
been studied quite extensively in the past, both experimentally and theoretically [[TJ[2]. Typically 
they display in the ground state and at low temperatures, the antiferromagnetic (AF) and, when 
increasing the field, the spin-flop (SF) phase, as depicted in Fig. 1. 

There may also occur intermediate or mixed spin structures of 'biconical' [ [3] (BC) type. 
They are characterized by distinct turning angles, @a and 6_b, with respect to the easy axis for 
neighboring spins situated on the two sublattices, see also Fig. 1. Quite a few experiments related 
to BC structures have been performed over the years [ HI [5] . 

Recent studies deal with the influence of defects [0(7], of the lattice dimension, of details of the 
anisotropic interactions, and of quantum effects on the phase diagrams. In particular, topological 
and (multi)critical properties of the phase diagrams have been elucidated. 

The aim of the present contribution is to draw attention to recent analyses on the prototypical 
XXZ Heisenberg antiferromagnet in a field and variants. The XXZ model is assumed to have a 



(a) (b) (c) 




Figure 1. Spin orientations on neighboring sites showing perfect antiferromagnetic (a), biconical 
(b), and spin-flop (c) structures. 
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uniaxial exchange anisotropy, with the magnetic field being along the easy axis. By then adding 
crystal-field single-ion terms of quadratic or cubic form, favoring non-uniaxial spin orderings, bi- 
conical structures may be stabilized. Much actual interest in these BC structures stems from the 
fact that they correspond to supersolid structures, when transcribing the Heisenbcrg antiferromag- 
net to a quantum lattice gas [ El [9j [TOl [Tl] . We shall consider both classical models, with spinvectors 
of length unity, and quantum models, with spin S=l/2 and S=l. In addition, to discuss the role 
of lattice dimension, results for chains, square lattices, and cubic lattices will be presented. 

2. XXZ antiferromagnets 

As a starting point of theoretical studies on uniaxially anisotropic Heisenberg antiferromagnets, 
one often considers the prototypical XXZ model, with the Hamiltonian 

-Hxxz = jJ2i A ^ S i + S ? S ^ + S i S l ] - H E S ' W 

i,j i 

where J > is the antiferromagnetic exchange coupling between spins being located on neighboring 
lattice sites i and j. A is the exchange anisotropy, 1 > A > in the case of uniaxiality, and H is 
the applied magnetic field along the easy axis, the z-axis. 

In the following, we shall discuss the phase diagrams of the XXZ antiferromagnets for square 
and cubic lattices, showing similar as well as quite distinct features. The phase diagrams have been 
obtained using mainly Monte Carlo (MC) simulations applying the Metropolis method for classical 
models and applying stochastic series expansion techniques for quantum systems, augmented by 
finite-size analyses. 

In our MC simulations we studied two-dimensional systems with up to 240 2 classical and up to 
150 2 quantum, S=l/2, spins. In three dimensions, we simulated systems with up to 32 3 classical 
spins. To estimate error bars, we averaged over several rather long runs. In all cases, full periodic 
boundary conditions have been employed. 

2.1. Square lattice 
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Figure 2. Phase diagram of the classical XXZ model on a square lattice near the maximum of 
the boundary line of the AF phase, A = 4/5, as obtained using Monte Carlo simulations. 



The classical XXZ model on a square lattice is known to display in the (temperature T, field 
H)-plane ordered AF and SF phases [[T2l [T3l [14l [T5] , see Fig. 2. The transition to the paramagnetic 
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phase belongs to the Ising universality class along the boundary line of the AF phase, and to the 
XY (or Kosterlitz-Thouless) universality class along the boundary line of the SF phase. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 3. Probability P2{Oa, Ob) showing the correlations between the tilt angles Oa and Ob on 
neighboring sites for the XXZ antiferromagnet on a square lattice, with 80 x 80 spins, at H/ J = 
2.41, UbT/J — 0.255, and A = g. pi is proportional to the gray scale. The superimposed black 
line depicts the relation between the two angles Oa and 0s in the ground state, eq. (2). 

Extensive Monte Carlo simulations augmented by finite-size analyses [ [13j [14] [15] suggest 
that there is no direct transition between the AF and SF phases, albeit the two ordered phases 
approach each other very closely at low temperatures, see Fig. 2. The intervening, extremely narrow 
paramagnetic (P) phase and the two ordered phases arise from the highly degenerate ground state 
occurring at the field H c = AJ\/1 — A 2 . At that point, not only AF and SF configurations have 
the same energy, but also biconical structures. The BC structures, Fig. 1(b), may be described by 
the tilt angles Qa and b , being uniquely interrelated by [ [15] 

9b = arccos ■ (2) 

\1 - VI- A 2 cose A J 

where the tilt angle Qa ranges from to tt. Indeed, the degenerate BC configurations have been 
argued to lead to the intervening disordered phase between the AF and SF phases [ [15] . Their 
presence at low temperatures may be demonstrated conveniently by monitoring the probability 
P2(Qa,Qb) to find the two tilt angles at neighboring sites, as illustrated in Fig. 3. Similar obser- 
vations hold for the anisotropic XY antiferromagnet on a square lattice [ [16] . 

The degenerate point at (T = 0, H c ) may be coined a 'hidden tetracritical point' [ [TBI IT?] . 

Due to the degeneracy, various physical quantities take at T = and H c nontrivial values, 
depending continuously on A. For instance, the magnetization per site is approximately 0.22 for 
A = 0.8 [UJ]. 

In contrast, the quantum XXZ antiferromagnet with S=l/2 on a square lattice seems to exhibit 
a direct transition of first order between the AF and SF phases [[TS]. The triple point, at which 
the AF, SF and P phases join, has been argued to be a critical endpoint, with the boundary line 
between the AF and P phases being of first order close to that triple point [ [TS] . The evidence 
had been provided mainly by quantum Monte Carlo simulations. We did quantum Monte Carlo 
simulations as well, for the same anisotropy, namely A = 2/3, using stochastic series expansions 
with directed loop updates [ [19], enlargening the system sizes and improving the statistics of 
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Figure 4. Positions of the maxima of the magnetization histograms as a function of the inverse 
system size, 1/L, for the S=l/2 XXZ antiferromagnet on a square lattice, with A = 2/3. The 
inset exemplifies two histograms for systems of linear size L = 32 (circles) and L = 150 (squares) 
at k B T/J = 0.13 and the coexistence fields H/J = 1.23075 and H/J = 1.232245 [ITS]. 



the simulations in comparison to the previous Monte Carlo study [ [T^l HH] ■ Our results suggest 
that the location of that special point may be shifted towards lower temperatures, compared to 
the previous estimate [[15]. This is illustrated in Fig. 4, where very accurate MC data for various 
system sizes are shown at a temperature closely above that of the proposed critical cndpoint, where 
the previous study suggested a transition of first order between the AF and P phases. Monitoring 
the size dependence of the peak positions in the magnetization histograms, Fig. 4, the two peaks, 
may well coincide in the thermodynamic limit, signaling a continuous transition [[H]. Because we 
do not find any evidence for BC structures in our simulations of the quantum model, there seems 
to be no mechanism for destroying long-range order at low temperatures in the quantum case. 
Thence the direct transition between the AF and SF phases seems to be allowed, in accordance 
with recent quantum MC findings [|5D]. In general, quantum fluctuations may substantially reduce 
BC structures, compared to the classical case, as will be discussed below. 

2.2. Cubic lattice 

The classical XXZ antiferromagnet on the cubic lattice displays also ordered AF and SF phases, 
see Fig. 5. In agreement with previous simulations [HI], we find the transition between the AF 
and SF phases to be of first order [ H2] . There is again a highly degenerate ground state, now 
at H c = 6JV1 — A 2 , with BC configurations as described by eq. (2). Indeed, by monitoring the 
probability p2 {<9a , 6s ) , BC configurations are observed to contribute to the thermal flucuations at 
low temperatures. But, close to the transition between the AF and SF phases, p2 shows pronounced 
local maxima corresponding to these two phases, indicating phase coexistence at a transition of 
first order. In addition, analysis of other standard thermodynamic quantities give evidence for a 
transition of first order [I22j. 

The nature of the triple point, at which the AF, SF, and P phases meet, has been subject to 
recent renormalization group (RNG) analyses, yielding either some kind of critical endpoint [ [23] 
or, possibly depending on the degree of anisotropy A, a bicritical point [[24]. The latter variant 
had been suggested in early RNG studies [|3] as well. 

From our MC simulations [H3 we can locate the triple point accurately at ksT / J — 1.025 ± 
0.015 and H/J = 3.90 ± 0.03, improving the previous estimate [[21]. We observe no evidence 
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for a first-order transition at the SF-P or AF-P boundary lines in the vicinity of the triple 
point at A = 0.8 [[22], being compatibel with the existence of a bicritical point. In fact, we find 
critical exponents, for instance for staggered longitudinal (i.e. along the easy axis) and transversal 
susceptibilities, at the AF-P boundary line to belong to the Ising universality class, while those 
at the SF-P boundary line are clearly consistent with the XY universality class. Perhaps, Monte 
Carlo studies at different values of the exchange anisotropy A may provide further insights to 
clarify the predictions of the RNG treatments. 



3. Adding single-ion anisotropics 

There is no thermally stable biconical phase in the classical XXZ antiferromagnet. BC config- 
urations occur in the ground state only at the special field, H c , separating AF and SF structures. 
BC fluctuations seem to lead to a narrow disordered phase between the AF and SF phases in two 
dimensions. In three (and, presumably, higher) dimensions they do not destroy the phase transition 
of first order between the AF and SF phases, as predicted by mean-field theory. 

However, BC structures may be stabilized over a wide range of fields, when introducing, for 
example, in the XXZ antiferromagnet additional quadratic or cubic single-ion anisotropics favoring 
non-uniaxial spin orientations. In the following we shall consider first phase diagrams, in the (T, 
H) plane, for classical Hcisenbcrg antiferromagnets with quadratic anisotropics on square lattices 
and with cubic anisotropics on cubic lattices. Finally, we shall present results on the ground state 
phase diagram of a quantum, 5 = 1, Heisenberg chain with uniaxial exchange anisotropy plus a 
quadratic single-ion anisotropy in a field. 

We use straightforward ground state considerations and Monte Carlo techniques for the classical 
models. The quantum spin chain at zero temperature is analysed by applying density matrix 
renormalization group (DMRG) methods. Whenever feasible, the finite-size behavior is taken into 
account. 

In our MC simulations we study systems with up to 240 2 spins in two dimensions, and with up 
to 32 2 spins in three dimensions. Details on the ground state DMRG calculations for the 5=1 
quantum Heisenberg chain will be given below. 

3.1. Classical antiferromagnets in two and three dimensions 

In the case of the square lattice, we add to the XXZ model, eq. (1), a single-ion term of the 
quadratic form 
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H D =DY,(SD 2 (3) 

i 

where the single-ion anisotropy may, depending on the sign of D, enhance the uniaxial exchange 
anisotropy A (0 < A < 1), when D < 0, or it may introduce a competing planar anisotropy, D > 0. 
The full Hamiltonian is then H = Hxxz + T~iv- 
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D/J 

Figure 6. Ground state phase diagram of the classical Heisenberg model with exchange, A, and 
quadratice single-ion, D, anisotropy on a square lattice, for A = 0.8. The filled circle denotes 
the highly degenerate point at D — and H = 2.4J. 

From a straighforward analysis of the ground states one may obtain the phase diagram at zero 
temperature [ HH1 US] . An example is shown in Fig. 6. There the exchange anisotropy has been 
fixed to be A= 0.8. The resulting ground state structures are of AF, SF or BC type. Of course, at 
sufficiently large fields, all spins will eventually point in the direction of the field, i.e. one encounters 
the ferromagnetic spin configuration, giving rise to the P phase at non-zero temperatures. As for 
the XXZ antiferromagnet, the xy-components of the spins order antiferromagnetically in the SF 
and BC configurations, as depicted in Fig. 1, having rotational symmetry. 

The most interesting feature of the ground state phase diagram is the region of stable BC 
structures, evolving, at D > 0, from the highly degenerate point of the XXZ model. In that region, 
the degeneracy is lifted, and the tilt angles 6^ and 0#, being still uniquely interrelated, change 
continuously with the field H [[TBI Above a critical positive value of D, D/J > 2 — 2A, at 
vanishing field, the spins are oriented perpendicular to the easy axis. Applying then a field, a SF 
structure will be stable, as shown in Fig. 6. 

At D > 0, the BC ground states give rise to a stable biconical phase in between the AF and 
SF phases, as exemplified in Fig. 7, setting D/J = 0.2. In the BC phase, the order parameters of 
the AF and SF phases, the staggerd longitudinal and transversal magnetizations, do not vanish 
[H M EH [23 . The BC-AF transition is expected to belong to the XY (or Kosterlitz-Thouless) 
universality class, while at the BC-SF transition the order parameter of the AF phase vanishes in 
an Ising-like manner. 

In the BC phase, the dominant tilt angles Qa and Qb change continuously with the field, when 
fixing the temperature. That behavior may be inferred from the ground state properties. It has 
been confirmed in simulations by monitoring the probability p(Q) of encountering a spin with the 
tilt angle O at an arbitrary lattice site, showing two pronounced peaks corresponding to the values 
of Qa and 9s at given temperature and field in the BC phase [ [25] . 
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Figure 7. Phase diagram of the XXZ antiferromagnet on a square lattice with a competing 
single-ion anisotropy, A = 0.8 and D/J — 0.2. 



It seems an open problem whether the AF, SF, BC, and P phases meet at one point, presumably 
a tetracritical point [[55]. At present, we can locate it, e.g. for the case depicted in Fig. 7, only 
with a rather large uncertainty [ [TJd], [55] . 

At D < 0, BC structures are squeezed out by the single-ion term in the ground state. At low 
temperatures, there seems to be a direct transition of first order between the AF and SF phases [ 

na (la- 
in the case of the cubic lattice, we add to the XXZ model, eq. (1), a cubic single-ion term [ [55] 

w CA = F^[(sr) 4 + (^) 4 + (^) 4 ] (4) 

i 

where F denotes the strength of the cubic anisotropy. The sign of F determines whether the spins 
tend to align along the cubic axes, for F < 0, case 1, or, for F > 0, case 2, in the diagonal 
directions of the lattice. Because of these tendencies, the BC, (i.e. BC1 or BC2), structures, as 
well as the SF structures, show no full rotational invariance in the xy-plane, in contrast to the 
XXZ antiferromagnet, with or without quadratic single-ion term. Now, as sketched in Fig. 8, the 
discretized spin projections in the xy-pl&ne favor four directions [[22]. 



(a) (b) 




Figure 8. Sketch of discretized biconical structures for the Heisenberg antiferromagnet with cubic 
anisotropy, in the case of (a) F > 0, BC2, and (b) F < 0, BC1. 



The resulting ground state phase diagram of the full Hamiltonian, TL = 7ixxz + Wca, with 
fixed exchange anisotropy, A = 0.8, and varying cubic term, F, may be determined numerically 
without difficulty [[55 ] 157 ] . as depicted in Fig. 9. 

For F < 0, the transitions to the BC1 structures are typically of first order, with a jump in the 
tilt angles, 0, with respect to the z-axis, characterizing the BC configurations. However, in the 
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Figure 9. Ground states of the three-dimensional classical XXZ antiferromagnet with exchange 
anisotropy, A = 0.8, and varying the cubic term F. 



reentrance region between the SF and BC1 structures at F/J close to -1, the change in the tilt 
angles seems to be smooth [[22]. 

Obviously, at non-zero temperatures, several interesting scenarios leading, possibly, to multi- 
critical behavior, where AF, SF, BC, and P phases meet, may exist. So far, we focussed attention 
on two cases [ [22] : 

(i) Positive cubic anisotropy F > 0, at constant field H/ J = 1.8 [ 22J, see Figs. 9 and 10. At small 
values of F, there is an AF ordering at low temperatures. Above a critical value, F c — 0.218. ..J, 
the low-temperature phase is of BC2 type, followed by the AF and P phases, when increasing 
the temperature. The transition between the AF and P phases is found to belong to the Ising 
universality class, while the transition between the BC2 and AF phases seems to belong to the XY 
universality class, with the cubic term being then an irrelevant perturbation [ 1221 128j . 
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Figure 10. Phase diagram of the XXZ antiferromagnet, A = 0.8, on a cubic lattice with a 
positive cubic anisotropy F at fixed field, H/J =1.8. 



(ii) Negative F, fixing the cubic term, F/J = —2, and varying the field, see Fig. 9. In accordance 
with the ground state analysis, we observe, at sufficiently low temperatures, AF, SF, BC1, and 
P phases when increasing the field. Interestingly, the BC1 phase seems to become unstable when 
raising the temperature, with the other phases being still present [[52]. This may suggest that the 
three boundary lines between the BC1-P, SF-BC1, and SF-P phases meet at a multicritical point. 
Of course, further clarification and a search for other, possibly multicritical scenarios at different 
strengths of the cubic term, F, are encouraged. 
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3.2. Ground state phase diagrams of S=l anisotropic Heisenberg chains 

In the following, we shall consider the spin-1 anisotropic antiferromagnetic Heisenberg chain 
with a quadratic single-ion term in a field described by the Hamiltonian 

n = Y^(J(sfSf +l + svs? +1 + As?sf +1 ) + D{stf - bsd (5) 

i 

where i denotes the lattice sites. The antiferromagnetic exchange coupling J is again positive. 
For A > 1, there is an uniaxial exchange anisotropy, along the direction of the field, B > 0, the 
z-axis. Note the slight change in defining the exchange anisotropy A compared to that for the 
XXZ model as given in eq. (1). We here follow the standard notation of previous studies on closely 
related quantum spin chains [ [29l [301 131] ■ For the same reason, the field is now denoted by B. 
Depending on the sign of D, the single-ion term leads to a competing planar anisotropy, D > 0, 
or to an enhancing uniaxial anisotropy, D < 0, as stated already above. We shall analyse ground 
state properties using DMRG techniques [ [32J, [33] for chains with open boundary conditions and 
up to L = 128 sites. In addition, we determine the ground states of corresponding infinite chains 
with classical spin vectors of length one [ O [55] . 

We shall compare our findings on the quantum spin chain to ground state properties of the 
corresponding classical chains, eq. (5). To characterise the quantum spin structures, the analogue 
of the classical spin-flop configuration will be called "spin-liquid' (SL) structure, and the analogue 
of the classical biconical configuration will be called supersolid (SS) structure, following the actual 
standard terminology. 

We first briefly deal with the case D/J = A/2 [ [30l [3T] . The ground state phase diagram in 
the (A,B/J) plane has been found [[30] [31] to comprise AF, SL, SS, ferromagnetic (F), and (10), 
with a magnetization plateau at half saturation, phases. At small values of A and small fields, the 
Haldane phase is observed [ [30] [34] . In comparison, the corresponding classical spin chain shows 
a much broader biconical (BC) phase, being effectively replaced not only by the supersolid phase 
but also, largely, by the SL and (10) phases. The classical (10) phase becomes stable only in the 
limit of an Ising antiferromagnetic chain with a single-ion term, the Blume-Capel model [ [31] . 




Figure 11. Magnetization profiles in the spin-liquid phase for A = 2D/J = 3.5 at various total 
magnetizations for the 5* = 1 chain, eq. (5), with L — 63 sites. 



Interesting information is given by the magnetization profiles, mi — (5f), with brackets, (...), 
denoting quantum mechanical expectation values. In particular, we observe distinct profiles in the 
SS phase in between the AF and (10) and in between the AF and SL phases, respectively. For odd 
L, in the SS phase on approach to the (10) phase, the local magnetizations rrn at odd sites stay 
close to one, while at even sites they tend roughly to zero. In contrast, in the SS phase on approach 
to the SL phase, the magnetizations on odd and even sites tend to take on the same values [131]. 
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In the SL phase we find also (studying situations with A exceeding w 2.5) two distinct types 
of profiles, for finite chains, when varying M/L, where M is the total magnetization [[3T]: For 
M < L/2, the profiles exhibit a broad plateau in the center of the chain, as expected for a classical 
spin-flop (SF) structure, while pronounced modulations in rrii occur at M > L/2, see Fig. 11. 
This may signal a change from commensurate (C) to incommensurate (IC) structures [ 31j, with 
an interesting quantum phase transition. The suggestion is confirmed and quantified by analysing 
the Fourier transform of the profiles, especially at A = 3.5. The modulation in the IC region of 
the SL phase is described nicely by the wavenumber (setting the lattice spacing equal to one) 
q = 2ir(l — m), m = M/L. Such an IC modulation, with algebraic decay, is expected to hold in 
the entire SL phase of the spin-1/2 anisotropic Heisenberg chain in a field [[35]. In the classical 
variant, we find no IC structures in the SF phase, for finite and infinite chains. 
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Figure 12. Ground state phase diagram of the classical infinite anisotropic Heisenberg chain, eq. 
(5), with A = 5.0. 

Let us now turn to the case of fixed exchange anisotropy, A = 5.0, varying the single-ion term, 
D. The ground state phase diagrams for the classical and the spin— 1 chains are depicted in Figs. 
12 and 13, using DMRG calculations for chains with up to 63 sites for the quantum case, and 
straightforward ground state considerations [ [25] (checked by Monte Carlo data) for the infinite 
classical chain. We considered positive and negative single-ion anisotropies, D, finding, especially, 
several intrigiung quantum phase transitions. 
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of D/J — A/2, the broad BC phase of the classical chain is effectively replaced, in the quantum 
chain, by the corresponding, rather narrow SS phase as well as SL and (10) phases. 

In contrast to the case D/J = A/2, the supersolid phase is always bordered by the AF and SL 
phases. Accordingly, we observe only one type of magnetization profile. Illustrative examples are 
depicted in Fig. 14, at D / J = 3.0 and various fields. In the SS phase, at given single-ion anisotropy, 
D, and field, B, the magnetization takes on different values at odd and even sites in the center of 
the chain. The local magnetization m, tends to acquire the same value at odd and even sites on 
approach to the SL phase in the quantum chain. Actually, as described above, the classical BC 
phase is characterized by two tilt angles Qa and Qb for the two sublattices formed by neighboring 
sites, with the tilt angles approaching each other when getting closer to the SF phase [ IT5l I25j. 
Obviously, this behavior is completely analogous to the one depicted in Fig. 14. 
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Figure 14. Magnetization profiles in the supersolid phase at A = 5.0, D/J = 3.0, and fields B/J 
(corresponding to total magnetizations M) between 6.5 and 6.8. Spin chains with L — 63 sites 
are considered. The lower (upper) symbols denote the even (odd) sites. 



Increasing D, the 'large-D' phase [ [37J (Ml GUI] may eventually be stable, leading to other 
intrigiung quantum phase transitions. It corresponds to the planar phase in the classical model 
with vanishing field, with the spin vectors pointing perpendicular to the z-axis, being the ground 
state for D/J> 4. The new phase may be expected to give rise to a SL phase at non-zero fields. 
A discussion of this interesting aspect is, however, beyond the scope of the present contribution. 

For D < 0, biconical and supersolid phase are no longer stable. In fact, only AF, SL (or SF), and 
F structures are encountered, in accordance with the discussion in the preceeding subsection. The 
SL phase can be either commensurate, with a wide plateau in the magnetization profile away from 
the boundaries, or incommensurate, with modulations superimposed on the average magnetization. 
Obviously, there are two IC phases, see Fig. 13. The one, denoted by TC' in that figure, occurring 
essentially in between the AF and F phases, has been found before, having exponentially decaying 
transversal spin-spin correlations [ [29] , in contrast to the usual spin-liquids with algebraic decay. 
It has no analog in the classical model, see Fig. 12. The related transition between the IC and 
SLc phases has been obtained before, being either of first order or, at large average magnetization, 
continuous [[29]. Our results agree with that description. We find another C-IC transition line 
between the SLc and SLjc phases at somewhat larger fields for given D/J, see Fig. 12 (full 
circles). Note that this transition seems to take place at M/L significantly larger than 1/2 for 
D < 0. Increasing D, D > 0, in the vicinity of the (10) phase, the line goes over to the above 
discussed scenario with commensurate, M/L < 1/2, and incommensurate, M/L > 1/2, structures. 
By further increasing the planar single-ion anisotropy, D/J being larger than roughly 3.6, we 
observe in the 'SLc' phase close to the (10) phase modulated structures for rather short chains, 
L < 31. The possible, additional C-IC border is not displayed in Fig. 13. Indeed, a more detailed 
analysis of that region, taking into account finite-size effects, is desirable. 
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4. Summary 

In this contribution, we have presented results of recent and current studies on a variety of 
anisotropic Heisenberg antiferromagncts in a field, elucidating the role of exchange and crystal 
field single-ion anisotropies, of lattice dimension, and of quantum effects. 

The prototypical classical XXZ Hcicnbcrg antifcrromagnet with uniaxial exchange anisotropy in 
a field along the easy axis is shown to exhibit highly degenerate biconical ground states at the field 
separating AF and SF configurations. The BC structures give rise to a narrow disordered phase 
intervening between the AF and SF phases at low temperatures for square lattices. In contrast, 
for cubic lattices the BC fluctuations do not destroy the direct transition of first order between 
the two ordered phases. Quantum effects seem to suppress the BC structures in the S=l/2 XXZ 
model on a square lattice, leading, presumably, to a transition of first ortder between the AF and 
SF phases. 

Thermally stable BC phases may occur when adding to the exchange anisotropy of the XXZ 
antiferromagnet, for instance, a quadratic or cubic single-ion anisotropy favoring a non-uniaxial 
ordering of the spins. Resulting phase diagrams for classical two- and three-dimensional magnets 
exhibit intriguing multicritical points. 

Studying ground state properties of S — 1 anisotropic Heisenberg chains with an additional 
quadratic single-ion anisotropy, we observe, for finite chains, two distinct types of spin-liquid 
as well as two distinct types of supersolid structures. The region of stability of the supcrsolid 
structures is substantially reduced, compared to that of the corresponding biconical structures in 
the classical case. Of course, there are no thermal phase transitions in one-dimensional magnets 
with short-range interactions. However, there are interesting quantum phase transitions at zero 
temperature. 
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